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abstract 
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for  peek  quasi-stitic  pressures,  *l°*do**  “V  * Cfiited  states  end  Europe  have 

blest  loading  within  e“olo*ttr**’  thea  displayed  graphically  according  to  r els- 
been  collected  and  analysed,  end  th«  display^  ,ruphs  presented, 

tionehips  derived  fro®  the  quesi-statio  preasure  versus 

along  with  appropriate  curve  fits, the  reduced  duration  versus 
the  ratio  of  charge  weight  to  encl  *|l#  iapul#0  v.r.u.  the  reduced 

the  reduced  pressure,  end  the  reauo 

preasure.  \ 

^ nomenclature 


■ surfaoe  area  of  enclosure 

- unit  of  force 

- unit  of  length 
. unit  of  time 

- enclosure  volume 

- unit  of  energy 

■ charge  weight 

■ sound  speed  of  air 

- exponential  decay  constant 
m functional  relationship 

_ functional  relationship 
m functional  relationship 

- specific  (gee)  impulse 
m pressure 
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tjjjj  * duration 

aeff  B r*tio  of  effective  vent  area  to  tot'l  enoloaure  surface  area 
Y * ratio  of  specific  heats 
n = nondimensional  teru 
v * nondimensional  duration 
0 " standard  deviation 

Superscripts ; 

- (bar)  ■ nondimensional  term 
a (carrot)  - calculated  quantity  from  curve  fit 


introduction 


The  loading  from  an  explosive  charge  detonated  within  a vented  or  un- 
vested structure  consists  of  two  almost  distinct  phases.  The  first  phase  is 
that  of  reflected  blast  loading.  It  consists  of  the  initial  high  pressure, 
short  duration  reflected  wave,  plus  perhaps  several  later  reflected  pulses  ar- 
riving at  times  closely  approximated  by  twioe  the  average  time  of  arrival  at 
the  chamber  walls.  These  later  pulses  are  usually  attenuated  in  amplitude  be- 
cause of  irreversible  thermodynamle  process,  and  they  may  be  very  oomplex  in 
waveform  becauae  of  the  complexity  of  the  reflection  process  within  the  struc- 
ture, whether  vented  or  unvented.  Maxima  fox  the  initial  internal  blast  loads 
on  a structure  can  be  estimated  from  scaled  blast  data  or  theoretical  analyses 
of  normal  blast  wave  reflection  from  a rigid  wall.  Following  the  initial  and 
secondary  shook  wave  reflections  from  the  internal  walls,  the  pressure  settles 
to  a slowly  decaying  level  — the  shook  wave  phase  of  the  loading  is  over. 

The  second  phase  of  a slowly  deoaying  pressure  is  a function  of  the  volume  and 
vent  area  of  the  structure,  and  the  nature  and  energy  release  of  the  explo- 
sion. 


The  process  of  reflection  and  pressure  buildup  in  either  unvented  or 
poorly  vented  structures  has  been  recognized  for  some  time,  dating  from  World 
War  II  research  on  effects  of  bombs  and  explosives  detonated  within  enclo- 
sures. However,  very  little  data  were  available  from  WWII  and  no  attempt  was 
made  to  understand  or  relate  the  physical  processes  until  1968  when  Weibull 
\ [1]  correlated  peak  quasi-static  pressure  versus  the  charge  weight  for  a ser- 

[ ies  of  experiments  with  TNT  detonated  within  a vented  enclosure.  More  recent- 

I>  ly,  study  of  these  pressures  hat  revived  beoause  of  interest  in  design  of 

vented  and  unvented  explosion  containment  chambers. 

i' 

A typical  time  history  of  pressure  at  the  wall  of  a vented  structure  is 
shown  in  Figure  1,  One  can  see  that  the  maximum  quasi-static  pressure  is 
quite  difficult  to  define  beoause  it  is  obsoured  by  the  initial  shook  and 
first  few  reflected  shooks.  Obviously,  several  reflections  must  ooour  before 
? irreversible  processes  attenuate  the  shooks  and  convert  their  energy  to  quasi- 

statio  pressure.  Therefore,  it  seems  inappropriate  to  oall  point  A in  Figure 
1 the  peak  quasi-static  prossure,  although  this  is  the  point  used  by  Kingery, 

| et  al.  [2]  to  oompare  with  code  predictions  from  Proctor  and  Filler  [3]  and 
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Figure  1.  Typical  Pressure-Tine  History  of  an  Explosion  in  a Tented  Structure 


the  Kinney  and  Sewell  equation  [4].  A better  approach  ia  to  allow  tone  time 
for  establiahing  the  maximum  pressure,  such  as  point  B in  Figure  1.  The 
smaller  the  vent  area  to  the  total  wall  area,  the  closer  the  pressure  at  point 
B will  be  to  that  at  point  A,  and  in  the  limit  of  the  vent  area  being  zero 
(that  is,  an  unvented  enclosure),  the  pressures  at  points  A and  B will  coin- 
cide. 


Figure  1 also  illustrates  another  problem  inherent  in  reduction  of 
vented  pressure  data:  accurate  determination  of  duration  of  this  pressure. 

When  the  pressure  approaches  ambitfnt,  the  shock  reflections  have  largely  de- 
cayed. But,  the  pressure  approaches  the  baseline  nearly  asymptotically  so 
that  the  duration  is  quite  difficult  to  determine  accurately.  A possible  du- 
ration tg|tx  is  shown  in  the  figure. 

In  spite  of  complexities  in  the  venting  proceso,  gas  pressures  and 
their  durations  can  be  predicted  with  reasonable  accuracy,  particularly  if  one 
differentiates  between  these  relatively  long  term  and  low  amplitude  pressures 
from  the  internal  blast  pressures  resulting  from  blast  wave  impingement  and  re- 
flection. Generally  of  great  interest  in  the  blast  loading  of  structures  re- 
sulting from  interior  explosions  are  the  peak  quasi-statio  pressure,  the  dura- 
tion, and  the  total  impulse.  The  following  paragraphs  will  develop  and  de- 
scribe relationships  among  various  physical  parameters.  Functional  relation- 
ships will  be  derived  from  similitude  analysis,  but  the  exact  functional  form 
cannot  be  ascertained  from  this  analysis  without  invoking  restrictive,  simpli- 
fying assumptions.  However,  a large  quantity  of  experimental  data,  from  a va- 
riety of  sources  [1-22],  allow  empirical  relationships  to  be  obtained. 


SIMILITUDE  ANALYSIS 


A model  analysis  was  performed  to  determine  the  functional  form  of  the 
quasi-statio  pressure  versus  physical  parameters  pertinent  to  a vented  struc- 
ture, The  problem  is  envisioned  as  an  instantaneous  energy  release  of  magni- 
tude W inside  a confined  volume  V.  A vent  area  (ueff  A)  exists  through  which 
internal  gases  can  escape,  where  aeff  is  the  effective  ratio  of  vent  area  to 
total  oross-sectional  area  of  the  walls.  Ambient  atmospherio  pressure  p0  ex- 
ists initially  inside  and  outside  the  oonfined  volume.  To  define  an  oquation- 
of-state  for  the  gases  in  this  problem  requires  two  additional  parameters 
which  are  the  ratio  of  specific  heats  y and  speed  of  sound  a0.  Table  1 summa- 
rizes the  parameters  in  this  problem  and  lists  their  fundamental  dimensions  in 
a system  of  units  of  force  F,  length  L,  and  time  T. 

Nondimensional  numbers,  or  pi  terms,  oan  be  developed  from  the  list  of 
variables  iu  Table  1.  The  assumptions  in  this  analysis  are  all  in  the  defini- 
tion of  the  problem.  Phenomena  are  not  considered  which  do  not  have  parame- 
ters listed  in  the  table.  Probably  the  major  assumption  invoked  is  that  ther- 
mal effects  are  ignored  — in  other  words,  the  pressures  dissipate  through 
venting  and  not  through  the  conduction  of  heat  into  the  walls  of  the  struc- 
ture. An  acceptable  set  of  pi  terms  which  result  are: 
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Note,  however,  that  ng  adds  no  new  information  since  the  impulse  can  be  ex- 
plicitly obtained  by  integrating  the  pressure  with  respeot  to  time. 

In  general  terms,  dimensional  analysis  states  that  the  functional  for- 
mat for  the  reduced  pressure,  «i,  is  given  by: 
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If  we  are  only  interested  in  predicting  the  peek  quasi-static  pressure,  the  re- 
sult will  not  depend  upon  time,  hence  the  functionsl  form  must  be  invariant 
with  respect  to  the  last  pi  term.  Likewise,  for  y a constant  (as  it  would  be 
for  air),  the  functional  form  will  not  depend  upon  y,  hence: 


s 


iir 

!’■ 


if 


P 


<o 


slL 


,2/3 


(8) 


where  j!  is  the  ratio  of  the  maximum  absolute  quasi-static  pressure  to  ambient 
pressure,  i.e.: 


and  pqs  is  the  conventional  gage  quasi-static  pressure.  Provided  the  flow 
through  the  vent  is  snail  relative  to  the  rate  of  energy  release,  the  maximum 
pressure  will  occur  before  significant  venting  has  transpired.  And  since  the 
ambient  pressure  is  essentially  an  invariant,  Equation  (8)  can  then  be  written 
for  the  maximum  quasi-static  pressure: 
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(10) 


The  blow-down  time,  or  duration,  can  also  be  expressed  as  a functional 
relationship  with  respect  to  the  other  pi  terms: 


Til.  Atf  a> 
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(11) 


But  it  just  has  been  shown,  if  the  maximum  quaai-statio  pressure  is  reached 
before  significant  venting  occurs,  that  the  last  term  1/ (p0V)  is  a function  of 
the  first  term,  p/p0.  And,  since  y is  an  invariant,  Equation  (11)  becomes: 
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Based  on  a theoretical  analysis  of  ohambor  venting  by  Owozarek  [22], 
Baker  and  Oldham  [24]  showed  that 
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In  physical  terms,  Expression  (14)  states  that  the  blowd  <wn  time  is  directly 
proportional  to  the  chamber  volume  divided  by  the  effec  ive  venting  area,  not 
an  unexpected  result.  Expression  (13)  thus  allows  us  to  simplify  Equation 
(12)  by  defining  a new  scaled  time  T: 
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Thus,  the  seeled  duration  Is  also  only  a function  of  the  reduced  pressure. 

JThe  last  relationship  to  obtain  is  a nondiaensional,  or  reduced,  im- 
pulse ij.  Figure  2 shows  a simplified  form  for  gas  vesting  pressures.  In 
this  simplified  form,  the  gas  venting  pressure  is  assumed  to  follow  the  solid 
curve  which  rises  linearly  from  time  zero  until  it  reaches,  at  time  tj,  a 
ourve  which  is  deoaying  exponentially  from  an  initial  maximum  value  of  p^, 
where  pi  Is  the  absolute  quasi-statio  pressure  at  time  t - 0.  The  decay  then 
follows  the  time  history 


p(t) 


P1  e 
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(17) 


until  it  reaches  ambient  pressure  p0  at  time  t » tmax.  The  exponential  decay 
is  shown  to  agree  well  with  experiment  (Kingery,  et  al.  [2],  and  Schumaoher, 
et  al,  [ 51 ) . The  oross-hatohed  area  under  the  overpressure  ourve  is  defined 
as  the  gat  impulse  ig,  and  is  found  by  integrating  Equation  (17)  with  respeot 
to  time : 
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(18) 


The  duration,  tatx,  will  be  obtained  from  Equations  (15)  and  (16).  Likewise, 
as  will  be  shown  later,  the  exponential  decay  factor,  o,  oan  be  written  in 
terms  of  tnxx  and  $ . As  stated  earlier,  the  impulse  oan  thus  be  found  ex- 
plicitly from  the  other  nondimenaional  relationships,  but  it  is  still  useful 
to  display  the  impulse  graphically,  Sinoe  the  impulse  depends  upon  the  dura- 
tion anu  the  pressure,  and  the  soaling  faotor  for  time  is  given  by  Expression 
(14),  then  a suitable  ohoioe  of  parameters  to  scale  the  impulse  is: 
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Figur.  2.  Simplified  Pre.sure-Timo  History  for  Venting 
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But  the  maximum  pressure  and  the  sealed  duration  are  functions  only  of  the  re- 
duced pressure,  hence 


1 

8 


= h 


(20) 


Equations  (10),  (16),  and  (20)  express  the  functional  relationships  be- 
tween the  various  physical  parameters.  The  next  section  will  empirically 
determine  the  functional  forms  via  curve-fitting  of  experimental  data.  How- 
ever, a brief  discussion  of  the  effective  vent  area  ratio,  aeff,  is  in  order. 


Venting  can  be  geometrically  quite  complicated  for  some  structures, 
particularly  those  structures  referred  to  as  suppressive  structures  whioh  of- 
ten have  three  to  six  wall  layers  with  various  staggered  venting  patterns  so 
fragments  will  not  esoape  the  confinement.  For  multi-walled  confinement,  an 
affective  vent  area  ratio  must  be  computed.  To  compute  aeff  for  a multi- 
walled  structure,  we  have  used 


eff  i 


N i 

E h 


(21) 


where  N is  the  number  of  elements  in  a suppressive  structure  panel.  Although 
no  proof  of  this  relationship  is  presently  possible,  it  does  reaoh  the  appro- 
priate limits  for  small  and  large  numbers  of  plates.  For  example,  if  only  one  s 

plate  is  present,  a#fj>  - as  it  should.  If  an  infinite  number  of  plates  is 

present,  uef f •*  0,  witn  the  flow  completely  choked.  If  one  of  the  plates  is 
solid,  and  thus  has  a zero  a,  a0j f ■ 0 as  it  should.  If  all  plates  have  the 
same  value  for  a,  atf{  * a/N,  which  is  a number  smaller  than  a for  a single 
plate  as  would  be  expeoted.  In  each  member,  a is  defined  according  to: 


a 


A 

vent 

Awall 


(22) 


For  plates,  the  meaning  of  this  definition  is  obvious.  However,  in  angles  and 
louvres,  the  definition  is  less  obvious  since  snglea  and  louvres  are  more 
efficient  in  constricting  flow  than  are  plates  vith  holes.  Details  for  comput- 
ing the  a's  for  more  complicated  geometries  can  be  found  in  Baker  and  Oldham 
[24]  and  Baker,  et  al.  [25] . 


GRAPHICAL  DISPLAY  OF  DATA 


The  proceeding  discussion  determined  which  physical  parameters  are  in- 
terrelated. This  analysis  permits  the  Judicious  choice  of  parameters  to 
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display  experimental  data  graphically.  Equation  (10)  states  that  the  peak 
quasi-static  pressure  is  only  a function  of  the  charge  ener gy-to-volume  ratio. 
Thus,  a plot  of  experimental  data  of  pQg  versus  W/V  will  determine  this  func- 
tional relationship.  Figures  3 and  4 display  the  data  from  177  tests,  and  as 
oan  be  seen,  the  experimental  data  range  over  several  orders  of  magnitude. 
Figure  3 is  a graph  of  the  data  in  metric  units,  while  Figure  4 is  the  identi- 
cal graph  except  it  is  displayed  in  English  units.  The  data  include  tests 
conducted  with  the  following  high  explosives  (HE):  TNT,  PE1N,  PBX-9502,  50/50 

Pentolite,  dynamite,  C-4,  Comp  B,  and  RDX. 

One  approximation  has  been  made  in  plotting  these  data.  For  any  given 
high  explosive,  the  oharge  energy  is  directly  proportional  to  the  explosive 
mass.  Also,  the  energy-to-mass  ratio  for  most  high  explosives  is  approximate- 
ly the  same.  Figures  3 and  4,  for  convenience,  use  the  mass  of  the  explosive 
for  the  symbol  W.  No  attempt  has  been  made  to  normalize  all  the  high  explo- 
sives to  TNT  since  the  scatter  in  data  from  experiments  with  the  same  high  ex- 
plosive often  masks  the  effects  of  slight  variations  of  the  energy-to-maaa  ra- 
tio differences  between  explosives.  (For  oarefully  controlled  experiments, 
the  differences  in  effects  of  energy  variation  between  explosives  can  be  mea- 
sured. Indeed,  for  a aeries  of  experiments  oonduoted  by  Hokanson,  et  al. 

[21],  where  quasi-static  pressures  were  measured  for  bare  explosive  and  the 
same  explosive  mass  encased  in  plastic  and  aluminum,  the  contribution  of  the 
oxidized  oasing  to  the  peak  quasi-statio  pressure  was  theoretically  oomputed 
and  measured  experimentally.)  It  should  be  noted,  however,  that  if  explosion 
scenarios  other  than  HE  detonations  are  of  interest,  e.g.,  various  fuel/air 
mixtures,  then  the  TNT  equivalent  weight  should  first  be  determined  and  then 
used  for  W if  these  graphs  are  to  be  used  to  determine  the  peak  quasi-statio 
pressure. 

It  is  reasonable  to  expect  the  peak  quasi-statio  pressure  to  be  direct- 
ly proportional  to  the  charge  weight,  and  examination  of  the  data  in  Figure  3 
confirms  this  supposition  for  small  and  large  W/V.  For  intermediate  values  of 
W/V,  a transition  region  is  evident.  For  W/V  < 0.4  kg/m^,  complete  oxidation 
of  the  explosive  occurs.  But  if  W/V  is  too  large,  insufficient  oxygen  is 
available  to  convert  all  the  potential  energy  available  in  the  explosive 
oharge,  and  the  energy  release  is  reduced  by  the  ratio  of  the  heat  of  detona- 
tion to  the  heat  of  combustion.  Thus,  for  W/V  > 11.0  kg/m^,  the  primary  oxi- 
dizer available  is  that  in  tho  explosive  itself.  A transition  region,  0.4  < 
W/V  < 11.0,  connects  the  two  regions.  Linear  least-squares  curve  fits  have 
been  performed  on  the  data  in  the  two  end  regions,  and  are  shown  in  Figure  3.* 
A seventh-order  polynomial  of  the  form: 


* All  curve  fittings  in  this  article  have  been  performed  in  log-log  coordi- 
nates, Linear  thur,  refers  to  the  form  of  the  curve  in  a log-log  plot. 


log  P«c  ■ * + s (W/V)  + a,  [log  (W/V)]2  + a,  [log  (W/V)]3  + 

US  o 1 2 i 

t4  [log  (W/V)]4  + »s  [log  (W/V)]5  + t6  [log  (W/V)]6  + 

•7  [log  (W/V)]7  (23) 


was  then  used  to  ourve  fit  the  entire  range  of  data,  where  log  repreaenta  the 
logarithm  to  the  baae  10.  Sueh  a high  order  polynomial  ia  required  beoauas  of 
the  number  of  oonatraiata:  the  alopea  and  interoepta  at  each  end,  the  pointa 

where  the  polynomial  connecta  with  the  atraight  linea,  and  the  oonatraiat  that 
the  ourve  be  a leaat-aquarea  fit,  Thia  ia  a total  of  seven  oonatraiata  which 
atipulatea  at  loaat  a sixth-order  polynomial  for  the  ourve  fit.  But  the  ap- 
pearance of  the  data  in  Figure  3 implies  that  an  odd  function  (as  opposed  to 
an  even  function)  should  be  used.  Hence,  a seventh-order  polynomial  becomes 
the  minimum  order  polynomial  stipulated.  Table  2 lists  the  coefficients  of 
Equation  (23)  as  well  as  the  linear  expressions  for  the  two  asymptotes.  Table 
3 lists  the  comparable  coefficients  and  asymptotes  for  pressure  in  psi  and  W/V 
in  lb/ft^  (Figure  4). 

The  standard  deviation  for  Equation  (23)  has  also  been  computed  but 
needs  to  be  interpreted  properly.  The  standard  deviation,  o,  is  usually  used 
as  an  estimate  of  the  scatter  in  data,  or  error  in  predictions.  One  standard 
deviation  enoompassea  approximately  68  percent  of  all  data  values.  The  uncer- 
tainty in  estimating  an  observable  is  often  written  as  the  calculated  quantity 
plus  or  minus  one  standard  deviation: 


log  p 


QS 


lo*  pqs  i 
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(24) 


where  pq$  it  the  estimated  quasi-static  pressure,  pQg  is  the  computed  quasi- 
static  pressure  from  the  curve  fit,  and  oc  is  the  computed  standard  deviation. 
Define  o0  such  that 


log  oo  - ag 

The  right-hand  side  of  Equation  (24)  can  then  be  written  as 

log  pQS  + log  oo 

so  that  Equation  (24)  can  be  rewritten  as: 


(25) 
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, 3 

W/V  i 0.4  kg/a 

$QS  * 2.347  (Yl/V)0,8395 

a “1 *143 
0 

3 

W/V  l 11.0  kg/m 

$Qg  “ 1.1004  (W/V)0, 9202 

a m 1*300 
0 
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Table  3. 


Summary  of  Peak  pQg  Veraua  W/V 
[pal  Veraua  lb/ft^] 


log  Pne  “ 3.3138  + 0.952133  log  (W/V)  + -0.023074  [log  (W/V))  + 

QS 


0.317807  [log  (W/V))3  + 0.149364  [log  (W/V))4  + 0.374595  [In  (W/V)]5  + 
0.161978  [log  (W/V)]6  + 0.021486  [log  (W/V)]7 


Correlation  Coefficient,  r:  0.993 


One  Standard  Deviation:  o0  ■ 1.247 


n iftS,,.  j _ / i 247  n 

1.247  * PQS  1 X*  47  PQS 


Asymptote* : 


W/V  i 0.025  lb/ft 

J 


pQS  - 3495.  (W/V); 


0.8435 


a - 1.143 
o 


W/V  l 0.70  lb/ft 


QS 


2049.  (W/V) 


0.9393 


a - 1.300 
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A 
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Tablet  2 and  3 give  a0,  aa  well  at  the  correlation  coefficient  (whioh  it  a 
measure  of  the  confidence  of  the  curve  fit). 

Here  we  wonld  like  to  mention  our  uncertainty  of  whether  the  slopes  of 
the  two  asymptotes  should  have  the  same  value.  The  slopes  are  not  appreciably 
different,  and  particularly  with  the  scatter  prevalent  for  large  K/V,  it  is 
not  unreasonable  to  speculate  that  the  slope  should  be  identical.  However, 
for  the  present,  we  have  elected  to  report  the  linear  least-aquares,  that  is, 
tho  best  fit  to  the  data. 

In  deriving  Equation  (10),  the  assumption  was  made  that  the  flow 
through  the  vent  (for  acff  4 0)  is  small  relative  to  the  rate  of  energy  re- 
lease so  that  the  maximum  quasi-static  pressure  occurs  before  significant 
venting  has  transpired.  Keenan  and  Tancreto  [22]  obtained  no  measurable 
quasi-static  pressure  for  values  of  (aeff  A/V^^)  £ 0.772.  Of  the  data  plot- 
ted in  Figures  3 and  4,  the  maximum  reduoed  vent  area  ratio  was  0.3246.  Thus, 
Figures  3 or  4,  and  Tables  2 or  3,  are  valid  for 


i 0.3246 


(29) 


For  a vented  enclosure  (oeff  4 0),  Equation  (16)  suggests  that  the  du- 
ration, be  plotted  versus  the  reduoed  pressure,  p,  given  by  Equation  (9). 
Seventy  of  the  data  points  from  Figure  3 or  4 represent  vented  enclosures,  and 
these  are  plotted  in  Figure  3.  It  oan  be  seen  that  the  duration  has  consider- 
able scatter  beoause  of  the  difficulty  in  determining  when  the  overpressure 
hss  returned  to  ambient.  Note,  however,  that  the  uncertainty  in  duration  has 
negligible  impact  on  the  impulse  since  the  total  area  under  the  pressure  time 
curve  is  not  sensitive  to  the  exact  location  of  tm#x.  A linear  least-squares 
curve  fit  has  been  performed  on  the  dsta  and  is  shown  in  Figure  5,  and  the  re- 
sults are  summarized  in  Table  4. 


From  the  linear  ourve  fit,  it  is  straightforward  to  compute  tmtx  from 
Equation  (15)  as  a function  of  p: 


max 


a a ..  A 
o ef  f 


(0.4284) 


!0S Ho 


.0.3638 


(30) 
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The  constant  c in  Equation  (17)  can  now  be  evaluated: 


P 


o 


"C 

<»QS  + ® 


t 

BAX 


(31) 


0 “ 


(32) 


The  specific  impulse  it  thon  obtained  from  Equation  (18),  which  after  some  re- 
arrangement of  terms,  reduces  to: 


i 

> 


(M. 

\ln  p 


(33) 


where  p is  given  by  Equation  (9)  and  tBax  is  given  by  Equation  <30). 

At  we  have  already  stated,  and  ju'-t  nhown  with  Equation  (33),  the  spo- 
oifio  impulse  can  be  obtained  direotly  from  the  peak  quasi-static  pressure  and 
the  duration.  However,  because  of  the  interest  in  specific  impulse  for  com- 
puting the  loading  of  structures,  it  is  often  convenient  to  have  a graphical 
representation  of  specific  impulse.  Equation  (20)  indicates  that  an  appropri- 
ate parameter  for  the  abscissa  is  the  reduced  pressure.  Sufficient  informa- 
tion was  reported  to  compute  specific  impulses  for  73  of  the  tests.  Figure  6 
displays  these  reduced  impulses,  ig,  versus  the  reduoed  pressure.  A linear 
least-squares  curve  fit  was  also  performed  on  these  data,  and  is  displayed  in 
Figure  6.  Table  3 summarizes  the  curve  fitting  information. 

Quadratic  least-squares  curve  fits  were  also  performed  on  the  data  in 
Figures  S and  6.  However,  the  standard  deviations  differed  by  less  than  seven 
percent  between  the  linear  and  quadratic  curve  fits  for  duration,  and  differed 
by  only  two  percent  for  the  reduoed  impulso.  A two-sample  comparison  of  vari- 
ance was  performed  using  the  F ratio  test  at  s 99  percent  confidence  level. 

For  the  linear  and  quadratic  curve  fits  to  be  different  statistically,  the  ra- 
tio of  their  respective  o's  must  exceed  approximately  1.7.  Since  the  ratio  of 
their  o's  is  much  less  than  1.7,  there  is  no  significant  difference  in  the 
linear  versua  quadratic  curve  fits  — hence,  only  the  expressions  for  the  lin- 
ear curve  fits  have  been  reported. 


SUMMARY 


A sizeable  quantity  of  data  have  been  compiled  and  analyzed  to  obtain 
peak  quasi-static  pressure,  and  the  duration  and  impulse  for  explosions  within 
structures.  Similitude  analysis  indicated  an  appropriate  ohoioe  of  parameters 
for  graphically  displaying  the  data.  Peak  quasi-static  pressure  was  found  to 
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be  e function  of  charge  weight  to  chamber  volume.  Alto,  a nondimensional  dur- 
ation and  a nondimentional  specific  impulae  were  found  to  be  functions  only  of 
the  reduced,  i.e.,  nondimensional , pressure.  The  data  range  over  several  or- 
ders of  magnitude  and  have  thus  been  presented  on  log-log  plots.  Least- 
squares  curve  fits  have  been  performed  and  reported,  with  their  standard  devi- 
ations, to  provide  appropriate  analytic  functions  to  relate  the  physical  pa- 
rameters. Thus,  for  high-order  detonations  within  enclosures,  the  peak  quasi- 
static  pressure,  and  if  the  enclosure  is  vented,  the  duration  and  specific  im- 
pulse, can  be  ascsirtained  from  the  graphs  on  the  respective  analytic  expres- 
sions. 
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